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ABSTRACT 


This report covers the period of the grant from January 1991 until its ex- 
piration in June 1992. Together with an Interim Report [Ref. 9], it summarizes 
the research conducted under NASA Grant NAG9-357 on the topic “Applica- 
tion of Lanczos Vectors to Control Design of Flexible Structures.” The research 
concerns various ways to obtain reduced-order mathematical models of complex 
structures for use in dynamics analysis and in the design of control systems for 
these structures. This report summarizes researd^described in the following 
reports and papers that were written during the second half of the grant period. 
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• Su, Tzu-Jeng and Craig, Roy R., Jr., “Control Design Based on a Linear 
State Function Observer,” Ref. [3]. 
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rithm For Damped Structural Dynamics Systems,” Ref. [4]. 


ii 



Book Chapters 


• Su, Tzu-Jeng and Craig, Roy R., Jr., ‘Krylov Vector Methods for Model 
Reduction and Control of Flexible Structures,” Ref. [5]. 

Papers Presented at Technical Meetings 

• Su, Tzu-Jeng and Craig, Roy R., Jr., “Structural Control Design Based 
on Reduced-Order Observer,” Ref. [6]. 

• Su, Tzu-Jeng and Craig, Roy R., Jr., “An Unsymmetric Lanczos Algo- 
rithm For Damped Structural Dynamics Systems," Ref. [7]. 

• Craig, Roy R., Jr., Su, Tzu-Jeng, and Kim, Hyoung M., “Some Experi- 
ences with Krylov Vectors and Lanczos Vectors,” Ref. [8]. 

Bound Reports 

• Craig, Roy R., Jr. and Su, Tzu-Jeng, “Application of Lanczos Vectors to 
Control Design of Flexible Structures,” Ref. [9]. 

• Su, Tzu-Jeng and Craig, Roy R., Jr., “Unsymmetric Lanczos Model Re- 
duction and Linear State Function Observer for Flexible Structures,” 
Ref. [10]. 



TABLE OF CONTENTS 


ACKNOWLEDGEMENT i 

ABSTRACT « 

TABLE OF CONTENTS iv 

1. INTRODUCTION 1 

2. ABSTRACTS OF TECHNICAL PAPERS 3 

2.1 Krylov Model Reduction (Refs. 1, 2, 5, 10) 4 

2.2 Unsymmetric Lanczos Algorithms for Damped Structural Dy- 
namics Systems (Refs. 4, 7, 10, 12) 11 

2.3 Control Design Based on a Linear State Function Observer (Refs. 

3, 6, 10) 16 

2.4 Substructure-Based Control of Flexible Structures 22 

BIBLIOGRAPHY 28 


IV 



Chapter 1 


INTRODUCTION 


Together with Ref. 9, this report summarizes the research accomplished 
over a three-year period on the topic of the application of Krylov vectors and 
Lanczos vectors to the control of flexible structures. The goals stated for re- 
search under the grant are: 

1 . To develop a theory of reduced-order modeling of general linear systems 
based on the use of Lanczos vectors, and to apply the theory to the 
modeling of flexible structures. 

2. To address numerical issues that arise in the application of Lanczos vec- 
tors to reduced-order modeling, for example, sensitivity to choice of start- 
ing vectors, loss of orthogonality, etc. 

3. To develop control system design techniques employing Lanczos modeling 
of the controlled and residual systems, considering relevant issues such 
as stability of the closed-loop system, spillover, robustness, and compu- 
tational requirements. 

4. To apply Lanczos-based control system design to several typical prob- 
lems, for example, optimal co-located velocity feedback, dynamic output 
feedback, optimal control of finite-time slewing of a beam, etc. 
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Algorithms for reduced-order modeling of general linear systems and of 
damped and undamped structures (Goal 1) are described in Refs. [1,2,4,7,10,11, 
14, and 15]. Steps that may be taken to address several numerical issues that 
arise in the application of Lanczos vectors to reduced-order modeling (Goal 2) 
are described in Refs. [7,10-12]. References [1,3,6,10,11,13-15,20] address the 
topic of control system design techniques (Goals 3 and 4). Summary papers 
and reports covering research performed under the present grant appear in 
Refs. [5,8,9,16]. 

Chapter 2 contains abstracts of the major research publications and pre- 
sentations not previously summarized in the Interim Report [Ref. 9]. 



Chapter 2 


ABSTRACTS OF TECHNICAL PAPERS 

In Ref. [9] , which covers research accomplished during the first eighteen 
months of this grant, abstracts were presented on the following topics: 

• Controller Reduction by Preserving Impulse Response Energy 

• Substructuring Decomposition and Controller Synthesis 

• A Review of Model Reduction Methods for Structural Control Design 

• Recent Literature on Structural Modeling, Identification, and Analysis 

Abstracts of publications describing research conducted during the last eighteen 
months of the grant period are presented in this chapter. 


3 
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2.1 Krylov Model Reduction (Refs. 1, 2, 5, 10) 

Some studies have shown that Krylov/Lanczos-based reduced-order mod- 
els provide an alternative to normal-mode (eigenvector) reduced-order models 
in application to structural control problems. The formulation based on Krylov 
vectors can eliminate control spillover and observation spillover while leaving 
only the dynamic spillover terms to be considered. 

An undamped structural dynamics system can be described by the input- 
output equations 

Mx + Kx = Pu X&FC 1 , ueR 1 
y = Vx + Wx y € R m 

Model reduction of structural dynamics systems is usually based on the Rayleigh- 
Ritz method of selecting an n x r transformation matrix L such that 

x = Lx ( 2 ) 

where x € R r (r < n) is the reduced-order vector of (physical or generalized) 

coordinates. Then, the reduced system equation is given by 

M'x + Kx = Pu , . 

_ _. ( 3 ) 

y = Vx + Wx 

where M = L T ML , K = L T KL, P = L T P, V — V L, and W = W L. The 
projection matrix L can be chosen arbitrarily. However, by choosing L to be 
formed by a particular set of Krylov vectors, it can be shown that the resulting 
reduced-order model matches a set of parameters called low-frequency moments. 

For a general linear system 

z=Az + Bu zeR?,ueR l 

y = Cz y e R m 


( 4 ) 



5 


the low-frequency moments are defined by CA~ l B, i = 1,2,..., which are the 
coefficient matrices in the Taylor series expansion of the system transfer func- 
tion. Applying the Fourier transform to Eq. (la) yields the frequency response 
solution X(u) = (K - u 2 M)' 1 PU (u) , with X(u) and U{u) the Fourier trans- 
forms of x and u. If the system is assumed to have no rigid-body motion, then 
a Taylor expansion of the frequency response around u = 0 is possible. Thus, 

OO 

X(u) = (I - u 2 K~ 1 M)~ l K~ 1 PU(u) = £ u 2i (K~ 1 M) i K~ 1 PU (u>) (5) 

i— 0 

Combining Eq. (lb) and Eq. (5), the system output frequency response can be 
expressed as 

OO 

Y{u) = Y^lViK-'MYK-'P + juWiK-'MyK-'Pp^Uioj) (6) 

i=0 

In these expressions, V(K~ l M) l K~ 1 P and W \K~ l M) x K~ l P play roles similar 
to that of low-frequency moments in the first-order state-space formulation. To 
obtain the reduced-order model of Eqs. (3) let 

span {L} = span {Lp Ly L\y} (7) 


where 

L p =[k~ 1 P C K~ l M)K~ l P ••• (K~ l MyK~ l P ] 
L v =[k~ 1 V t (K- l M)K~ l V T ■■■ (K- 1 Af)*K- l V T ] (8) 
L w =[k~'W t (K~ 1 M)K- 1 W t ••• (K- 1 M) a K~ 1 W T ] 

for p, q, s> 0. Then the reduced system matches the low frequency moments 
V{K- l MyK- l P for i = 0, 1, . . . , p + q + 1 and W (K~ l M) { K~ l P , for i = 
0, 1, 2, ..., p + s+ 1. 
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The Lp matrix above is the generalized controllability matrix, and the 
Ly and Lw matrices are the generalized observability matrices of the dynamic 
system described by Eq. (1). The vectors contained in Lp are Krylov vectors 
that are generated in block form by 

Q\ = K~ l P 
Qi+\ = K-'MQi 

The first vector block, K~ l P, is the system’s static deflection due to the force 
distribution P. The vector block Q,+\ can be interpreted as the static deflection 
produced by the inertia force associated with the Qi. If only the dynamic 
response simulation is of interest, then L = Lp would be chosen. In this case, 
the reduced model matches p + 1 low-frequency moments. As to the vectors 
in Ly and Lw, a physical interpretation such as the “static deflection due to 
sensor distribution” may be inadequate. However, from an input-output point 
of view, Ly, Lw, and Lp are equally important as far as parameter-matching 
of the reduced-order model is concerned. 

Based on Eq. (7), the following algorithm may be used to generate a 
Krylov basis that produces a reduced-order model with the stated parameter- 
matching property. 

Krylov/Lanczos Algorithm 

(1) Starting block of vectors: 

(a) Q o = 0 

(b) Ro = K~ l P, P = linearly-independent portion of [P V T W T ] 

(c) KqKRq = UqHqUq ( singular-value decomposition ) 
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(d) Q i = RoUo'Eo 2 ( normalization ) 

(2) For j = 1, 2, . . . , fc — 1, repeat: 

(e) = K~ x MQj 

(f) R : = - QjVlj - Qj-\Bj (orthogonalization) 

Aj = QjKTl,, B j = 

(g) RjKRj = UjEjUj ( singular-value decomposition ) 

(h) Q j+ i — RjBj+i = RjUjT.J 5 ( normalization ) 

(3) Form the k-block projection matrix L = [ Qi Q 2 Qk ]• 

This algorithm is a Krylov algorithm, because the L matrix is generated by 
a Krylov recurrence formula (Step e). It is a Lanczos algorithm because the 
orthogonalization scheme is a three-term recursion scheme (Step f). 

If the projection matrix L generated by the above algorithm is em- 
ployed to perform model reduction, then the reduced-order model matches 
the low-frequency moments V {K~ l M)' K~ l P and W(K~ X M)' K~ l P , for i = 
0, 1, 2, ..., 2k — 1. It can also be shown that the reduced-order model 
approximates the lower natural frequencies of the full-order model. 

One interesting feature of the transformed system equation in Krylov 
coordinates is that it has a special form. Because of the special choice of start- 
ing vectors, K-orthogonalization, and three-term recurrence, the transformed 
system equation has a mass matrix in block-tridi agonal form, a stiffness matrix 
equal to the identity matrix, and force distribution and measurement distribu- 
tion matrices with nonzero elements only in the first block. The transformed 
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system equation has the form 


X X 


' X 1 

XXX 


0 

X X • 


o ; 

. . . 

X + X = i 

. 

• * X 


. 

X X 


0 


y=[x 0 0 ••• 0 ]x + [ x 0 0 ■■■ 0 ]i 
where x denotes the location of nonzero elements. This special form reflects the 
structure of a tandem system (Fig. 1), in which only subsystem is directly 



Figure 1. The Structure of a Tandem System. 


controlled and measured while the remaining subsystems, £*, i = 2, 3, ..., 
are excited through chained dynamic coupling. In control applications, this 
tandem structure of the dynamic equation eliminates the control spillover and 
the observation spillover, but there is still dynamic spillover. For dynamic 
response calculations, the block-tridi agonal form can lead to an efficient time- 
step solution and can save storage. 

The previous model-reduction strategy can be extended to damped struc- 
tural dynamics systems, which are described by the linear input-output equa- 
tions 
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Mx + Dx + Kx = Pu 
y = Vx + W x 

To arrive at an algorithm for constructing a reduced-order model that matches 
low-frequency moments, it is easier to start from the first-order formulation. 
The first-order differential equation equivalent to Eq. (10) can be expressed as 


D M 
M 0 

y=[V IV] 



( 11 ) 


or „ 

Mz + Kz = Pu 

y = Vz 

with 


( 12 ) 


M = 


D M 
M 0 


K 0 
0 - 



V = [ V w] (13) 


QU 1 1 


97+i J 




[ 9?1 


1 

© 

1 


(14) 


This leads to the following recurrence formula for the Krylov blocks: 

-K~ l D -K~ l M 
I 0 

Superscripts d and v denote displacement and velocity portions of the vector, 
respectively. The matrix containing the generated vector sequence is called a 
Krylov matrix. It has the form 


'of of of 
.05 Of Of 

Krylov subspaces that are generated by Eq. (14) and that have the above 
form produce a projection subspace L that has the desired moment-matching 
property. 
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' Qf Qt Qi ••• Qi 
. o Qi Qi ••• QJ_! 


(15a) 


be the sequence of vectors generated by Eq. (14) with K l P the starting block 
of vectors, i.e., Qi = K~ l P, Q\ = 0, and let 


Ly = 


pi pi pi 
m pf pi 


pi 

pi - 1 


(15b) 


be the subspace of vectors generated by Eq. (14) with K 1 V r the starting block 
of vectors, i.e., Pi = K~ 1 V T , P? = If the projection matrix L is 

chosen such that 


span {L} = span { Qi • • • Q d v Pi ■■■ Pi Pi } (16) 

then the reduced-order model of the damped structural dynamics system matches 
the system parameters V(k~ l M)'K~ 1 P, for i = 0, 1, . . . , p + q — 1. 

A 48-DOF plane truss structure is used in Reference [1] to demon- 
strate the Krylov model reduction method. It is shown that the Krylov-based 
reduced-models can approximate the system’s impulse response better than the 
normal-mode reduced-models. The same structure is also used to illustrate the 
efficacy of Krylov vector method in the application to flexible structure control 
problems. 
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2.2 Unsymmetric Lanczos Algorithms for Damped Struc- 
tural Dynamics Systems (Refs. 4, 7, 10, 12) 

The Krylov model reduction method described in the previous section has 
a restriction: the damping matrix has to be symmetric. Although most passive 
damping mechanisms yield a symmetric damping matrix, there are cases when 
the damping matrix is unsymmetric. For structures, unsymmetric damping 
may arise from either active feedback control or Coriolis forces. To deal with 
general, nonsymmetric damping, the usual approach is to write the system’s 
dynamic equation in first-order state-space form. Then, an unsymmetric Lanc- 
zos algorithm is used to create a basis for model reduction of the first-order 
differential equations. All the existing unsymmetric Lanczos algorithms are 
two-sided [12], In this section, a one-sided, three-term, block-Lanczos iteration 
scheme will be presented. This iteration scheme, together with a special choice 
of starting block of vectors, can be used to create a set of Lanczos vectors. The 
projection matrix formed by the set of Lanczos vectors transforms the system 
matrix A into an almost skew-symmetric, block- tridiagonal form. 

First, assume that there exists a set of blocks of vectors Q it i = 1, 2, . . . , k, 
that satisfy 

‘ ?i Gx 
Hi T* Qi 

A[Q X Q 2 ••• Qk] = [Qi Q2 Qt] H * \ (17) 

The above equation implies the following iteration formula 


AQi = Qi-iGi-\ + Qt^i + Qi+iHi 


(18) 
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Let us further assume that the Qi's are orthonormalized with respect to the 
inverse of the controllability grammian of the system. That is 

I if i = j 


QjW-'Qj = 


0 


if * 7 * 3 


(19) 


Then, by premultiplying Eq. (18) by Qj-\W C \ QjW c l , and Qf + l W c 1 respec- 
tively, and using the orthogonality conditions in Eq. (19), it can be shown that 


Gi—\ = QhW-'AQi 

Ti = QjW~ x AQi (20) 

Kt = QT+iW-'AQi 

Let the starting block of vectors of the iteration formula in Eq. (18) be the B 
matrix after being normalized with respect to W~ x . That is, let 


Q X = BU 0 T,^ with UoXoUf = B T W~ l B (21) 

where UqEoUq is the singular- value decomposition of B T W~ 1 B with UjUo = I 
and Eo being the diagonal matrix containing the singular values. With this 
special choice of starting block of vectors, the following two identities can be 
derived. 


T. + Fj = 


= -Gl 

(22) 

— E 0 if i = 1 

0 if i ^ 1 

(23) 


These two identities indicate that the block tridiagonal matrix in Eq. (17) is 
almost skew-symmetric, except the diagonal elements of T\. 


Let 


L= [ Q, Q 2 ■■■ Q„ } 


(24) 
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Tx Gi 

-Gf r<2 g 2 
-Gl ^3 


(25) 


Then, Eq. (17) can be written as 

AL = LT (26) 

The orthogonality conditions in Eq. (19) imply that L T W~ X L = I. Therefore, 
the % unsymmetric system matrix A can be transformed into an almost skew- 
symmetric, block-tridiagonal form as shown below 

A = L T W~ l AL = T (27) 

The unsymmetric Lanczos algorithm proposed in this section is described 
by the following algorithm. 


One-Sided Unsymmetric Lanczos Algorithm 
Given A, B, and C. 

( 1 ) Solve AW C + W C A T + BB T = 0 for W c . 

( 2 ) Starting vectors: 

(a) Qo = 0 

(b) Rq = B 

(c) RZW-'R* = UqLqUq = C/oaEoaC/oa (SVD) 

_ 1 

(d) Qx = RqUoqE 0 a ( normalization ) 


14 

(3) For i = l, 2, k — 1 repeat: 

(e) Ri = AQi - Qi-iGi-i - QiTi ( orthog .) 

Qi - 1 = QT-iW-'AQ, , 0 Fi = QfW-'AQi 

(f) RJW-'FU = UiZiUT = U ia Z iQ Ul ( SVD ) 

- 1 

(g) <3i+i = RiUi Q T, i( ? ( normalization ) 

end 

(4) Form the k-block transformation matrix L — [Q i Q 2 ■ Qk\- 

It is easy to show that the Lanczos vectors generated by the above algo- 
rithm have the following property 

span [Qi,Q 2 , - ,Qk] = span \Q\,AQ U ... ,v4* -1 Qi] 

The sequence of vectors in the right-hand side matrix of the above equation is 
frequently called a set of Krylov vectors, which can be generated by the simple 
iteration formula 

Qi+i = AQi 

The set of Lanczos vectors generated by the proposed algorithm is basically the 

orthonormal ized set of the Krylov vectors. Since the starting block of vectors 

_ 1 

is chosen to be Q\ = BU 0 T, 0 ^ , the Lanczos vectors created are in the column 
subspace of B, AB , A k ~ l B j, which, in linear system terminology, is 

called the controllability matrix of the system. Therefore, the projection matrix 
L defined in Eq. (24) has a column subspace that is the same as that of the 
system’s controllability matrix. 

The set of Lanczos vectors generated by the Lanczos iteration formula 
can be used as a basis for model reduction. Let the L matrix be partitioned 
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into 



where L R € /2 nxr corresponds to the retained portion and L T E j R nx ( n_r ) cor- 
responds to the truncated portion. Then, the state vector can be decomposed 
into 

z = L r z r T L t Zt 

The reduced system equation is described by 

z R = A r z r + B r u 
y — C R z R 

where 

A r = T r = L t r W- 1 AL r> b r = L T R W- l B , C R = cl r 

Only L r is required to produce system matrices of the reduced-order model. 

Since the matrices of the reduced system satisfy T R + T R + B R B R = 0, 
the controllability grammian of the reduced system is an identity matrix. Also, 
the reduced system in Eq. (29) is a completely controllable system because 
the Lanczos vectors are in the controllable subspace of the system. According 
to the Lyapunov theorem, if a controllable system has a positive definite con- 
trollability grammian, then the system is stable. Therefore, the reduced-order 
model obtained by the Lanczos model reduction method is a stable system. 

The advantages of the proposed one-sided, unsymmetric Lanczos algo- 
rithm over the other existing unsymmetric Lanczos algorithms are: (1) the 
numerical breakdown problem that usually occurs in applying the two-sided 
unsymmetric Lanczos method is not present, (2) the Lanczos vectors that are 
produced lie in the controllable and observable subspace, (3) the reduced-order 


(28) 

(29) 
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model is guaranteed to be stable, (4) a shifting scheme can be used for un- 
stable systems, (5) the flexibility of the choice of starting vector can yield 
more accurate reduced-order models, and (6) the method is derived for general 
multi-input/multi-output systems. A linear system example and a plane truss 
structure example are included in Ref. [4] to show the efficacy of the proposed 
method. 

Reference [12] describes a two-sided unsymmetric block Lanczos algo- 
rithm that generates a set of left Lanczos vectors and a set of right Lanczos 
vectors (analogous to sets of left eigenvectors and right eigenvectors). These 
two sets of Lanczos vectors form a basis that transforms the system equation 
to an unsymmetric block-tridiagonal form. Computational enhancements that 
are described in Ref. [12] produce a very robust two-sided algorithm. 

2.3 Control Design Based on a Linear State Function 
Observer (Refs. 3, 6, 10) 

A method to design low-order controllers for large-scale systems is de- 
rived from the theory of linear state function observers. The observer design 
problem is considered as the reconstruction of a linear function of the state 
vector. The linear state function to be reconstructed is the given control law. 
Then, based on the derivation for linear state function observers, the observer 
design is formulated as a parameter optimization problem. The optimization 
objective is to generate a matrix that is close to the given feedback gain matrix. 
Based on that matrix, the form of the observer and a new control law can be 
determined. 
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Consider a linear, time-invariant system of the form 

z = Az + Bu, z E R n , u E R l 
y = Cz, y E R m 

A linear state feedback control law has the form 

u — Kz 

where the l x n feedback gain matrix K is designed to achieve some spec- 
ified performance objective by using an existing control design approach, for 
instance, linear quadratic control theory or the eigenvalues/eigenvectors assign- 
ment method. Due to the limited output measurement, the full state feedback 
control law in Eq. (31) in general cannot be realized. In order to implement the 
state feedback control, an observer is required to reconstruct the complete state 
vector. It is well known that for an n-th order system with m outputs, it is 
always possible to construct a reduced-order observer of order n—m, with poles 
arbitrarily placed (subject to complex paring), to yield an asymptotic estimate 
of the states. But, for large scale dynamic systems with very few outputs, 
an observer of order n — m is usually too large for implementation purposes. 
Hence, it is desirable to design an observer with considerably reduced order r, 
r < (n — m). 

Assume that the feedback gain matrix K can be decomposed into the 

form 

K = GT , G E R lxr , T E R rxn (32) 

where l < r < n - m. Then, the feedback control law in Eq. (31) becomes 


(30) 


(31) 


u = GTz 


(33) 
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which can be implemented by using the measurement of the linear state function 
Tz. Now, consider a r-th order observer of the form 

q = Eq + TBu + Fy, qe R T (34) 

The observer state vector q is to approximate Tz in the following asymptotic 
sense 

lim [q(t) — Tz(t)] — 0 (35) 

It is apparent that the observer system matrices E and F cannot be arbitrary 
and that they must satisfy some condition(s). Define the error vector 

e — q — Tz (36) 

Then, from Eqs. (30) and (34), 

e = Eq + FCz - TAz (37) 

If E and F are chosen to satisfy the Lyapunov equation 

ET - TA + FC = 0 (38) 


then Eq. (37) becomes 

e = Ee (39) 

The error dynamics is governed by the stability of the E matrix. 

In summary, in order to realize the control law in Eq. (31) exactly, the 
observer described by Eq. (34) must satisfy the following three conditions: 


(1) E is a stable matrix, 

(2) ET - TA + FC = 0, and 
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(3) K = GT 


If a feedback gain matrix K is already designed and given, it is reasonable 
to choose E and F matrices in such a way that the T matrix obtained from 
the Lyapunov equation (38) is “close” to K. The closeness of two matrices can 
be defined by the following two-norm measure 


\\K_-cr [| 2 

7 II K lb 

If T is close to K, then there exists a G matrix such that 7 is small, 
matrix that minimizes 7 is 


(40) 
The G 


G= KT + 


(41) 


where T + is the Moore-Penrose pseudo-inverse of T. 

To determine E and F matrices so that they yield a T matrix as close 
to K as possible, a parameter optimization technique can be used. All of the 
elements of E and F are considered as optimization parameters. 

Although E € R rxr and F 6 R rxm , only r x (m + 1) parameters are 
needed to characterize the observer system. The form of the E matrix can be 
chosen to be block-diagonal 


E = diag (A*) 


(42) 


in which 

is the i-th block associated with the pair of complex conjugate poles — <7* ± 
ju>i. For first-order poles, the E matrix has negative real numbers on the 
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diagonal. With the E matrix chosen in the block-diagonal form, the Lyapunov 
equation (38) can be split into independent partitions as 


'A, 


" h ' 




'/i ' 

a 2 


1 2 

— 

t2 

A + 

h 


(44) 


or, 

A it{ — t{A + fiC = 0, i = 1,2,... (45) 


where U and /, are the z-th partitions of T and F corresponding to the z-th 
block of E. 


A block-by-block optimization procedure can be used to determine U and 
fi. For the z-th partition, the optimization problem is described as: 


Minimize 7(C i,<**,/i) = 


K - KT+Ti || 2 


where T, = 


Ti-i 

U 


II K || 3 

, with U being the solution of 


UA = A iU + ftC 
and subject to 


(46) 


Ctmin — Cl — Cl max J ^imin — ^ 


l max 


Starting out with T 0 = [0], the observer system matrices E and F are optimized 
block by block with the objective of obtaining a T matrix that is close to K. 
Although block-by-block optimization does not yield a global minimum, it has 
fewer optimization parameters and the computation cost is more economic than 
a global optimization approach. A general criterion for choosing the bounds 
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for £ and Ui is to place the observer poles such that the observation process is 
fast enough to provide estimated states for feedback. 

Based upon the above derivation, a semi-inverse observer-based control 
design procedure can be summarized by the following steps: 

(1) Determine a state feedback gain matrix K. Set i = 1, and T — [0]. 

(2) Set the bounds for Q and Ui. Solve the optimization problem (19) for U, 
Ai, and /*. 

(3) Expand the T matrix by adding in Check the stability of A + BGT 
{G = KT + ). 

(4) Determine if the closed-loop system is stable with acceptable perfor- 
mance. If yes, stop; else set i — i + 1, and go to Step (2). 

The above control design procedure is called a semi-inverse procedure 
because the final control law is not designed before the observer, but is a result 
obtained from the observer design. The observer design, however, is not com- 
pletely independent of the control design, because it tries to yield a control law 
that is close to a given design. The dimension of the observer does not have 
to be specified in advance. Hence, the proposed method offers considerable 
flexibility for the design of a linear state feedback control law and its associ- 
ated linear state function observer. As the T matrix expands in Step (3), the 
order of the observer grows. If, at certain point the stability property and the 
performance of the closed-loop system is satisfactory, the design procedure can 
stop there. The observer that is obtained has a dimension that is the same as 
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the number of row vectors in the T matrix. A four-disk system and a lightly- 
damped beam are used as examples in Ref. [3] to demonstrate the efficacy of 
the proposed control design method. 

2.4 Substructure-Based Control of Flexible Structures 

Earlier work on this topic is summarized in the Interim Report [Ref. 9]. 
Research on this topic has continued under the present grant and two compan- 
ion grants (See Refs. [17] and [18]). The topic of robustness of substructural 
controllers was studied in Ref. [17]. Reference 18 compares the SCS method 
of substructural controller synthesis, which was developed under the present 
grant, with the CCS method described by K. D. Young in Ref. [19]. 

This section describes a substructure-based controller design method 
called Augmented Physical Component Synthesis (APCS). This approach evolved 
from an earlier substructure-based method called Substructural Controller Syn- 
thesis (SCS) [21]. The SCS method is very versatile because it is based on the 
individual substructure models and never explicitly uses the complete structure 
model in the controller design; this way, it can readily acoomodate a changing 
structure. However, this same versatility can lead to a very suboptimal or un- 
stable global controller. The APCS method uses component modal models of 
all of the substructures comprising the complete structure to create augmented 
components which, while still spatially distinct, approximate the response of 
a specific part of the complete structure. For example, consider a structure 
composed of two substructures, a, /?. An augmented a component is created 
by appending a very low-order model of the adjacent /? component. This is 
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depicted pictorially in Fig. 2.1. 

a substructure (3 substructure 


Augmented a component 

Low order model of . 

a substructure 


Low order model of 
P substructure 


Augmented p component 
Figure 2.1: Augmented components. 


An augmented /3 component is created similarly. When there are more than two 
substructures, the augmented components are created by sequentially, coupling 
reduced order models of the other substructures and appending them to the 
substructure of interest. With this method, each component will experience 
all of the actuators and contribute to all of the sensors. This allows for the 
specification of a global performance objective, which is not possible with other 
substructure-based methods. Essentially, the method attempts to account for 
the influences of a substructure’s neighbors in the design of the substructure 
controller in order to reduce the suboptimality of the controller. 

All dynamic controllers that use measurement data consist of two parts: 
the regulator and the observer. These are usually designed independently. 
First, a state feedback gain matrix is obtained for each component by some 
method. For example, for a linear quadratic (LQ) controller the component 
regulator gain matrix is obtained by minimizing a performance index associated 
with the augmented component, i : 
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1 T rrt 

f = i / [** QV + u r Ku]<if (47) 

2 ./o 

Minimization of this performance index for each component produces the com- 
ponent’s optimal feedback gain matrix, G *. Matrix (7* can be partitioned as 


u = — [G a G a 



(48) 


to reflect the gains associated with the original substructure states, z l s , and 
with the appended states, z\. Equation (48) manifests the contribution of the 
augmented component’s states to all of the actuators on the complete structure. 
A global gain matrix will be created from the component gain matrices, so the 
appended states cannot participate in the feedback (i.e., G a = 0) and these 
gains are simply truncated. The component control is then given by 


= Of 


2 S 



(49) 


The states of each component contribute to every actuator, so consequently G i 
has no zero rows. Therefore, upon creating the global gain matrix, the gains 
associated with the appended states are replaced by the component gain ma- 
trices of the other components. The global gain matrix is created by combining 
the component gain matrices. For a two component example, the APCS gain 
matrix is obtained by 
G° 


G° 


G* ] ( T T Y 


(50) 


where the superscript I means a Moore-Penrose pseudoinverse. The use of 
(T T ) / instead of the substructure coupling matrix , T , has the effect of aver- 
aging the gains at the common interface rather that summing them. In the 
limiting case, when full-order models of the substructures are appended to 
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create augmented components, the APCS gain matrix G° is the optimal gain 
matrix because each augmented component is a full-order model. 

The global observer is designed analogously. The state estimator of a 

system represented by a state space model 

z — Az + Bu + Tw 

y = Cz + v 

is of the form 

i = Az + Bu + L{y - Cz) 
u = —Gz 

In the APCS method, an estimator gain matrix, L\ is designed for each com- 
ponent. Each matrix V can be partitioned according to 


( 51 ) 


( 52 ) 


V = 


L a 

La 


( 53 ) 


to indicate the effect of the error correction term, L{y — Cz ), on the estimates 
of the original substructure states and appended states. As in the regulator 
design, the gains associated with the appended generalized states must be 
removed before a global gain matrix can be formed. The global estimation 
gain matrix is formed by combining the component estimation gain matrices 


L° = T 1 


L a 

IP 


( 54 ) 


Again the pseudoinverse is used so that the gains at the interface states will be 
averaged. In the limiting case, the APCS estimator gain matrix is the same as 
the global optimal one. 

The complete structure controller for a linear time invariant structural 
dynamics system is 

z = (A - BG° - L°C)z + L°y 
u = —G°z 


( 55 ) 
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Matrices A, B, C are defined in Eq. (51) and matrices G° and L° are de- 
fined in Eqs. (50) and (54), respectively. Because no a priori method exisits 
to evaluate the closed-loop properties of the complete structure based on the 
closed-loop properties of substructures, stability is not guaranteed. However, 
since the augmented component open-loop properties approximate those of the 
complete structure and the component gain matrices are assembled, the quality 
of APCS controllers should be quite good. Results from a simulation study of a 
long spring-mass-damper chain are encouraging. The complete structure is cre- 
ated from three substructures of eleven nodes each as shown in Fig. 2.2. Some 
results are shown in Figs. 2.3a,b, where the performance of a global optimal 
LQG controller and the APCS controller are compared. The regulator objec- 
tive was to minimize the position and velocity at each node. The process noise 
and the measurement noise intensities were set at 1 x 10 -2 and 1 x 10 -4 , respec- 
tively. Simulation studies for a more representative example of a large flexible 
structrure are currently being conducted. The results of these simulations will 
be published in a future CAR report. 

12 10 u 

&v Gv vtwfi 

Substructure 

1 2 3 29 30 31 

Complete structure assembled from 3 substructures 
Figure 2.2: Structure used in APSC example. 
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Figure 2.3: Comparison of optimal and APCS controllers. 
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